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ABSTRACT 

Supermassive black hole dynamics during galaxy mergers is crucial in determining the rate of black hole mergers and cosmic black 
hole growth. As simulations achieve higher resolution, it becomes important to assess whether the black hole dynamics is influenced 
by the treatment of the interstellar medium in different simulation codes. We here compare simulations of black hole growth in galaxy 
mergers with two codes: the Smoothed Particle Hydrodynamics code gasoline, and the Adaptive Mesh Refinement code ramses. We 
seek to identify predictions of these models that are robust despite differences in hydrodynamic methods and implementations of 
sub-grid physics. We find that the general behavior is consistent between codes. Black hole accretion is minimal while the galaxies 
are well-separated (and even as they “fly-by” within 10 kpc at first pericenter). At late stages, when the galaxies pass within a few 
kpc, tidal torques drive nuclear gas inflow that triggers bursts of black hole accretion accompanied by star formation. We also note 
quantitative discrepancies that are model-dependent: our ramses simulations show less star formation and black hole growth, and a 
smoother gas distribution with larger clumps and filaments, than our gasoline simulations. We attribute these differences primarily 
to the sub-grid models for black hole fueling and feedback and gas thermodynamics. The main conclusion is that differences exist 
quantitatively between codes, and this should be kept in mind when making comparisons with observations. However, reassuringly, 
both codes capture the same dynamical behaviours in terms of triggering of black hole accretion, star formation, and black hole 
dynamics. 
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1. Introduction 

Galaxy mergers are thought to be transformational events in 
galaxy evolution. Mergers transform stellar disks into spheroids 
(e.g. |Toomre & Toomre||1972[ |Gerhard| 1981[ [Negroponte &| 
White! 1983[ [Barnes & Hemquist||1996| ). Via tidal torques, they 
tend to compress gas into the central regions of galaxies, trigger- 
ing powerful starbursts (e.g. [Sanders et al.|1988[[Bames & Hem-j 
quist||199r[ [Mihos & Hernquist!1996| ). The increase in nuclear 

gas is further thought to fuel growth in galaxies’ central super- 
ma ssive black holes, resulting in active galactic nuclei (AGNs , 
e.g. [Sanders et al.|1988t|Hemquist|1989] Di Matteo et al.|2005] ). 


[et al.|2005}[Rupke & Veilleux|2011t [Gabor & Bournaud|2014[ ). 

AGNs sometimes power radio jets that can heat intergalactic and 
intracluster gas (e.g. [Fabian et al.|200(H[Voit & Donahue|2005 


els of galaxy evolution as a primary actor in regulating the star 
formation ra tes (SFRs) and stellar masses of the most massive 
galaxies (e.g.[Hopkins et al. [20061 [Croton et al. [20061 [Somerville 
[etal.|2008| ). - 

Numerical hydrodynamic simulations have led to many ad¬ 
vances in our understanding of galaxy mergers and black hole 
fueling, but we have not yet developed a complete understanding 
of how results depend on details of these models. Different stud¬ 
ies have used various treatments of hydrodynamics and feedback 
processes (such as supernova and AGN output), as well as dif¬ 
ferent resolution. Details of the feedback treatment can have an 
important effect on black hole growth and its impact on galaxies 


Correlations between supermassive black hole mass and 
global galaxy properties (see [Kormendy & Ho[[2013[ and ref¬ 
erences therein) suggest a possible evolutionary link between 
black hole and galaxy growth (e.g. Silk & Rees[[1998[ [Wyithe[ 
[& Loeb[[2003[ ). The energetic output from AGNs may provide 
a physical driver for this link: in some cases, they emit suffi¬ 
cient energy to heat up all the cold gas in a galaxy, and they may 
trigger pow erful outflows (e.g. [Crenshaw et al.|2003[[Di Matteo[ 


(e.g.[Debuhr et al.|201 l[[Wurster & Thacker|2013b|al [Newton & 
[Kay|2Q13[ ). 

Differences in hydrodynamic method can also lead to differ¬ 
ences in simulated galaxy (and inter-galactic) properties in var¬ 
ious contexts. Recent work shows that Smoothed Particle Hy¬ 
drodynamics (SPH; [Lucy[[1977t [Gingold & Monaghan[[1977| ), 
historically a commonly-used method in extragalactic astro¬ 
physics for solving the equations of hydrodynamics, is inaccu¬ 
rate in certain circumstances - for example in resolving shocks 
and Kelvin-Helmholtz and Rayleigh-Taylor instabilities ( [Agertz| 
et al.[2007]). Grid technique s, such as Adaptive Mesh Refinement 


Randall et al. 201 1[ ). AGNs have been invoked in many mod- (AMR; Berger & Colella[|l989[ ), generally improve upon these 
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problems ( [Agertz et al.|2007| ), but with drawbacks including ad- 
vection errors, angular momentum conservation, and numerical 
diffusion (e.g. |Wadsley et al.|2008[|Hahn et al.|20Tm see discus¬ 
sion and references in Hopkins |2014|). More recen t hybrid tech¬ 
niques employing a moving mesh ( |Springel|2010] ) or no mesh at 
all ( |Hopkins|2014| ) help resolve these difficulties. 

In cosmological simulat ions, traditiona l SPH as imple¬ 
mented in the gadget code ( |Springel||2005| ) tends to allow a 
smaller quantity of gas to cool and fuel galaxies than mesh-based 
codes, leading to smaller gas disks (|Keres et aL]|2012t [Vogels- 1 


jberger et al.|2Q12t|Scannapieco et al.|2012| ). The apparently sup¬ 
pressed cooling in SPH relative to moving mesh methods also 
allows the formation of more prominent hot halos in idealized 
merger simulations, along with clumps and filaments in those 
halos ( [Hayward et al. |20ll) . This is true despite the fact that, in 
cosmological SPH simulations, a smaller proportion of gas ac¬ 
creting onto galaxies passes through a hot phase ( [Nelson et al. 


[2013[ ). (Note that recent improvements t o SPH have helped al¬ 
leviate many of these discrepancies; e.g. [Beck et al.[[MT5] and 
references therein. Such improvements have recently been im¬ 
plemented in GASOLINE as described in e.g. [Keller et~SI|[2Q14[ ) 
Despite these important differences, it appears that differences in 
feedback implementations cause more prominent changes than 
differences in the numerical method for hydrodynamics ( [Scan-| 
napieco et al.[[2Q12| ). Moving forward, large code-comparison 


projects like AGORA (Kim et al. 2014) should help elucidate 


2011). We use standard physical recipes in both codes for star 


these i ssues in the context of galaxy evolution. With the excep¬ 
tion of [Hayward et al.[(2014[ ), the impact of different codes and 
numerical techniques on supermassive black hole fueling and 
feedback has scarcely been studied. 

In this work, we compare high-resolution simulations of 
black hole growth in galaxy mergers with two codes: the gaso¬ 
line SPH code, and the ramses AMR code. While many authors 
have used SPH codes to study black holes in idealized galaxy 
mergers, relatively few have used AMR codes (see e.g. [Kim et al.[ 


formation, black hole growth, and stellar and black hole feed¬ 
back. In ^we describe the simulations. Then we highlight sim¬ 
ilarities and differences, and attempt to explain them, in ^ We 
conclude in ^ 

2. Simulations 

For our analysis, we focus on four representative simulations of 
galaxy mergers: two run with gasoline and two with ramses. The 
merging galaxies have mass ratios of 2:1 and 4:1, gas fractions 
of 30%, and disks oriented coplanar with their orbits. We have 
also run some resolution tests and simulations with other mass 
ratios and orbital configurations (see [Capelo et al.|2Q15[ ), which 
give similar results. Below we describe numerical details of the 
GASOLINE simulations - including details of the initial conditions 
- then those of the ramses simulations. 


2.1. Methods: GASOLINE 

We use a subset of the suite of gasoline merger simulations de¬ 
scribed fully in [Capelo et al. p015 ). We summarize the simula¬ 
tions here, gasoline ( Wadsley et al. 2004[ ), an N-body SPH code, 
is based on pkdgrav ( [Stadel|2Q01[ ), which uses a tree method to 
calculate gravitational dynamics among particles. We note that 
these simulations do not include recent improvements in the SPH 
method ( [Keller et al.|2014[ ). 

As a Lagrangian particle code, the resolution is automatically 
adaptive - the highest resolution occurs in the densest regions. 


Dark matter particles, initial star particles, and initial gas parti¬ 
cles have masses of 1.1x10^ M©, 3.3x10^ M©, and 4.6x10^ M© 
(respectively), with softening lengths of 30 pc, 10 pc, and 20 pc 
(respectively). 

GASOLINE includes models for gas cooling, star formation, su¬ 
pernovae and stellar winds, and black hole accretion and feed¬ 
back. The simulations use a standard model for gas cooling, 
which incorporates metal cooling but no heating by a UV ra¬ 
diation background ( [Shen et ^[2010[ ). A temperature floor of 
500 K is imposed. 

Star formation occurs as a random process (see [Katz|1992| ) 
in gas particles colder than 6000 K and denser than 100 particles 
cm“^. The SFR is calculated by assuming that a fraction 6* = 
1.5% of the eligible gas forms into stars per star formation time, 
where the star formation time is the greater of the gas free fall 
time or the gas cooling time. With a probability consistent with 
the star formation rate and timestep, a star-forming gas particle 
will convert some of its mass into a collisionless star particle. 

Stellar feedback includes a blast wave model for Type II su¬ 
pernovae ( [Stinson et al.|2006} . In this model, supernova energy 
is calculated based on the mass of young stars that should ex¬ 
plode, and this energy is distributed (as thermal energy) among 
gas particles that are neighbors of the star particle. To mimic an 
expanding blast wave, gas cooling is turned off for a time that 
depends on the local gas conditions and the feedback energy. 

Nuclear black holes are treated as sink particles that accrete 
surrounding gas ( [Bellovary et al.|20T0| ). The accretion rate onto 
the black hole is calculated separately for each neighbo r gas par¬ 
ticle based on a Bondi accretion rate (cf. [Bondi[[1952[ ), and the 
sum of the rates from each gas particle yields the total accretion 
rate. Mass is then removed from the gas particles in proportion 
to their contribution to the accretion rate, and the mass is added 
to the black hole. We note that the Bondi accretion rate her e in¬ 
cludes a boost factor of a = 3 (see [Booth & Schaye||2009[ for a 
discussion of boost factors). The accretion rate is capped at a 
times the Eddington Limit. A fixed fraction 6r = 0.1 of the ac¬ 
creted mass-energy is emitted as radiation at each timestep, and 
a coupling fraction 6/ = 0.001 of the radiated energy injected as 
thermal energy into the nearest gas particle. 


2.1.1. GASOLINE simulation set up and initial conditions 


We set up mergers of two disk galaxies as in Capelo et al. (2015). 
Galaxies begin with parabolic orbits ( [Benson|2005[ ), initial sep¬ 
arations the sum of the two galaxies’ virial radii, and separation 
at first pericenter equal to 20% of the virial radius of the larger 
galaxy ( Khochfar & Burkert[2006[ ). The galactic disks are copla¬ 
nar with their orbits, and both galaxies rotate in the same direc¬ 
tion as the orbit (prograde-prograde mergers). 

Each galaxy is modeled as a dark matter halo, stellar bulge, 
stellar disk, and gaseous disk ( [Springel & White] 1999[ [Springel 
et al.[2005] ), plus a supermassive black hole. The dark matter halo 
follows an NEW ( [Navarro et al.|199^ profile up to the virial ra¬ 
dius, with an exponential decay beyond. It has spin parameter 
A = 0.04 and concentration c = 3. Th e stellar bulge m akes up 
0.8% of the virial mass, and follows a [Hemqui^ ( [ 1 990\ profile. 
The galactic disk makes up 4% of the virial mass, and follows 
an exponential surface density profile. The disk scale radius is 
derived from conservation of angular momentum of the mate¬ 
rial making up the disk. Thirty percent of the disk mass is gas, 
so that the fraction of total baryons in gas (the common obser¬ 
vational definition of gas fraction) is 0.3(0.04Mvir)/(0.04Mvir + 
O.OOSMyir) = 0.25. Superm assive black holes with masses 2x 


10 ^ times the bulge mass ( Marconi & Hunt|[2003 ) are placed 
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Fig. 1. Snapshots showing gas surface density for our 4:1 ramses (top) and gasoline (bottom) simulations. Left panels show both galaxies at first 
apocenter, with arrows roughly indicating the direction of travel of the secondary galaxy. Middle panels show a zoom-in on the primary galaxy at 
first apocenter, and right panels show both galaxies at second pericenter. In the right panels, circles mark the positions of the black holes (in each 
simulation, the more massive black hole is to the left). 


at the centers of initialized galaxies. In each simulation, the pri¬ 
mary galaxy has a virial mass of 2.21 x 10^^ M©, bulge mass 
1.77 X 10^ Mq, disk mass 8.83 x 10^ M©, and disk scale radius 
of 1.13 kpc. The secondary galaxies have their masses scaled 
down by factors of two and four, respectively. 

Model galaxies undergo a relaxation period of 100 Myr in 
isolation before they begin the merger. During this period, the 
star formation efficiency 6>, is gradually increased to its final 
value, 1.5%, to prevent unphysical bursts of supernova feedback. 
After this relaxation, the supermassive black hole masses are re¬ 
set to their initial values, and the galaxies are placed in appropri¬ 
ate orbits for the mergers. 


2 . 2 . Mef/iods; RAMSES 


We have run a small new suite of merger simulations using the 
RAMSES AM R code (|Teyssier||2002|). Ou r set up borrows many 
aspects from [Gabor & Boumaud| ( MT3l|20I4| ) and jPerret et al. 
( 20I4| ), and mostly uses standard recipes for physical processes. 
RAMSES solves the equations of hydrodynamics on the mesh, 
while it treats collisionless matter (dark matter and stars) as par¬ 
ticles. The effects of gravity are calculated on the mesh using 
a multigrid method, and particle accelerations are interpolated 
using a cloud-in-cell method. 

We use the standard quasi-lagrangian mesh refinement crite¬ 
ria: a cell is refined if it contains > 30 dark matter particles, or 
if the mass in the cell exceeds 5 x 10^ M©. In addition, we re¬ 
fine cells whose mass is sufficiently large that the local gas Jeans 


length is not resolved by at least 4 cell widths ( jTruelove et al. 
\1991) . These refinement criteria ensure that the highest resolu¬ 
tion is applied to the densest regions, where star-forming clouds 
form. We allow refinement up to level 15, corresponding to a 
minimum cell size of about 7.6 pc (with a box side length of 
250 kpc). 


The thermodynamics is treated with a standard model for 
gas cooling, including metal cooling (e.g. [Teyssier et al.|20I3 ). 
We impose an overall temperature floor of 100 K, as well as a 
density-dependent temperature floor (the Jeans polytrope) that 
ensures the local Jeans length in the smallest grid cells is al¬ 
ways resolved by at least 4 cells ( jMachacek et al.||200I| ). The 
Jeans polytrope acts as an additional pressure (or pressure floor) 
that prevents numerical fragmentation. The normalization of this 
temperature floor depends on resolution, and is higher for lower 
resolution simulations. 


RAMSES includes models for star formation and supernova 
feedback. Star formation occurs as a random process in cells 
whose gas (particle number) density exceeds 100 H atoms cm“^. 
In cells above this threshold density, a fraction 6>, = I per cent of 
the gas is assumed to form into stars per free fall time, yielding 
an SFR for each cell. A new collisionless star particle is created 
in the cell with a probability based on the SFR and the timestep. 
New star particles inherit the position and velocity of the gas 
cell out of which they formed, but they are decoupled and will 
generally move into other cells. 

After a delay of 10 Myr, supernovae explode at the locations 
of newly formed star particles. Twenty per cent of the mass of 
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Fig. 2. BH separation in kpc (top), black hole accretion rate (BHAR) in M© yr“^ (second row), gas mass in M© within spheres up to 1 kpc 
in radius (third row), and SFR in M© yr“^ within the same spheres (bottom row) vs. time for a 4:1 merger with gasoline (left) and ramses 
(right). For each simulation, we show quantities for the primary galaxy/BH on the left, and the secondary galaxy/BH on the right. We show gas 
masses and SFRs wit hin rad ii of 200 pc to 1 kpc, in increments of 200 pc. Quantities are smoothed on 10 Myr timescales. Note that time-axes 
differ by ~ 20% (see g3.2.1| ). Vertical dotted lines mark local minima in the separation between the two black holes in all panels. The gasoline 
and RAMSES simulations show qualitatively similar results: enhanced accretion and star formation rates at second pericenter, and especially at 3rd 
pericenter/coalescence. 


the initial star particle is assumed to explode as supernovae, and 
for each 10 M© of exploding supernova, 10^^ ergs of thermal en¬ 
ergy is added to the hosting cell. Following [Stinson et al.| ( |20Q6] ) 
and jTeyssier et al.| ( |2Q13| ), cooling is delayed for 20 Myr in the 
supernova-heated cell to enable more efficient feedback. 


Supermassive black holes are represented by collisionless 
sink particles. The sink particles accrete gas according to a Bondi 
accretion rate (cf. |Bondi|1952| ), where the gas density and tem¬ 
perature are computed from a weighted averag e of all gas cells 
within 4Ax (where Ax is the smallest cell size, [Krumholz et al.j 
|2004| ). We use the standard Bondi formula, without a “boost” 
factor a ( [Booth & Schaye||2009| ). The accretion rate is capped 
at the Eddington limit, assuming a radiative efficiency 6^ = 0.1. 
Black hole particles merge once they pass within 4Ax of one an¬ 
other. 


Black hole feedback is implemented as a thermal deposition 
of energy ( [Dubois et al.[[2010[ [Teyssier et al.[[2011| ). A fraction 
6^ = 0.1 of the accreted mass is assumed to convert into ra¬ 
diative energy, and a fraction 6^ = 0.15 of the radiative energy 
is assumed to couple to the surrounding gas as thermal heat- 
ing ([Booth & Scha^[2009t [Dubois et al. [[2010[ [Teyssier et al. 


201 1| ). Following Booth & Schaye[ ( 2009| ), we only inject AGN 
feedback energy if it is sufficient to heat the surrounding gas to 
Tmin,AGN = 10^ K. This prevents the deposited energy from being 
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immediately radiatively cooled from dense gas, leading to more 
efficient feedback. On timesteps where the feedback energy is 
insufficient to heat the gas to rniin,AGN, we store the feedback en¬ 
ergy to be added to that during the following timesteps. This stor¬ 
age is repeated until enough energy is stored to reach rniin,AGN 5 at 
which time it is released. We deposit the feedback energy in gas 
cells within 4Ax with a weighting where colder, denser gas ac¬ 
quires more of the energy ( [Gabor & Boumaud|2013[ ). If the post¬ 
injection gas temperature would exceed Tmax,AGN — 5 X 10^ K, 
then we iteratively expand the injection radius by 25 per cent to 
dilute the injection energy and lower the injection temperature. 
This maximum injection temperature prevents extremely high 
temperatures whic h can cause computational problems ( [Gabor 
[& Boumaud|2Q13| ). 


As noted in [Gabor & Bournaud] ( [2013| ), numerous tests show 
that black holes sometimes scatter from the centres of their host 
galaxies by hundreds of pc. In merger simulations, the scattering 
is even more pronounced. Some of this scattering is a numerical 
artifact - it essentially does not occur in the gasoline simula¬ 
tions (except on small scales, which can be reduced with im¬ 
proved dynamical friction modeling; see [Tremmel et al.[[2015 ). 
The scattering in ramses is reduced, but not eliminated, when 
a more accurate direct-summation N-body gravity solver is ap¬ 
plied to the sink particle (cf. [Bleuler & Teyssier][2014[ ). On the 
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other hand, some scattering (at the ~ 10^ pc level) is phy sically 
realistic given the dynamics in the nucleus (see |3.2.4| for ad¬ 
ditional discussion). To limit scattering, we adopt a well-known 
solution of assigning a black hole dynamical ma ss that is much 
larger than the true mass (e.g. |Debuhr et al.|2Q10| ). For the gravity 
calculation, the black hole is assigned a mass of 10^ M©, while 
we use the true mass for black hole-specific physics (e.g. Bondi 
accretion). This effectively keeps the black holes in their galac¬ 
tic centers by increasing the dynamical friction due to stars and 
dark matter. 


2.2.1. RAMSES simulation set up and initial conditions 


Initial conditions for the ramses simulations are generated sepa¬ 
rately from those for the gasoline runs, following the description 
in [Gabor & Bourna^ ( |2013| ). We use the same masses, density 
profiles and numbers of particles (for stars and for dark matter) 
as in the gasoline simulations, and initialize the positions and ve¬ 


locity distribution s with the Cartesian grid code from jBoumaud 
|& Comb^ ( |2002| ). 


The initial conditions generator determines equilibrium posi¬ 
tions and phase-space velocity distributions for stellar and dark 
matter particles, taking into account the gas disk contribution 
to the gravitational potential, and the gas density distribution if 
initialized analytically in ramses. We also must initialize gas 
cells outside the disk: technically, null densities could not be 
handled properly. We do this by setting their density very low 
(~ 10“^ cm“^) so that the initial contribution to the circum- 
galactic halo is marginal: this halo will form from outflows dur¬ 
ing the simulation. Supermassive black holes are added to the 
galactic centers at the start of the ramses runs. 

A necessary difference with the gasoline initial conditions 
is that the dark matter halo must be truncated to fit within the 
chosen simulation box size. The truncation radius, depending 
on the halo and box sizes, is 70 to 85% of the virial radius 
/?200- To maintain the halo truncation over time, we truncated 
the Maxwellian distribution of velocities for dark matter parti¬ 
cles at 90% of the local escape velocity (computed at the initial 
position of each particle), otherwise the halo expands spatially 
over time and dynamically evaporates from the ramses simula¬ 
tion box. When truncating the haloes, we made the choice to 
keep the same total halo mass, so at very large separations during 
the early phases of the interaction the two galaxies have the same 
total masses and keep similar orbits with the two codes, and un¬ 
dergo their collision under nearly identical configurations. Yet 
at lower separations once the systems start to overlap the halo 
densities start to differ: haloes generated this way end up with 
slightly higher masses within ^200 and smaller radii than the 
GASOLINE haloes (by a few percent typically), so short-range grav¬ 
ity and dynamical fraction are somewhat stronger in the ramses 
simulations. This will be discussed later when comparing the de¬ 
tailed evolution of each merging system. 

As in the gasoline case, we wish for each galaxy to undergo a 
relaxation period of ~ 100 Myr. We do not have a simple method 
to insert a fully relaxed galaxy into a ramses merger simulation, 
so the relaxation must occur during the main simulation. Given 
our orbital parameters, the galaxies remain quasi-isolated and 
undisturbed by tidal forces during the first 200 - 300 Myr of the 
merger, which is sufficient to allow relaxation (1-2 disk dynam¬ 
ical times) and allow the formation of internal features such as 
bars and spiral arms before the interaction itself occurs. 

For the merger, galaxies begin with the same orbital parame¬ 
ters (initial separation and velocities) and orientations (prograde- 


prograde) as in the gasoline simulations. As noted above, the 
mergers evolve in a cubic box with a size of 250 kpc. 


3. Results 


3.1. General behaviour: common features 

Figure presents face-on images of gas surface density in our 
4:1 RAMSES and gasoline simulations. First we show first apocen- 
ter, then a zoom-in of the primary galaxy at the same snapshot, 
and finally the merging galaxies at second pericenter. In Figure 
we compare the time evolution of various quantities during a 4:1 
merger in our gasoline and ramses simulations. At the top, we 
show the separation between the primary and secondary black 
holes. Below the separation plots, separately for the primary and 
secondary, we show the black hole accretion rates (BHARs), gas 
masses within 1 kpc of the black hole, and SFRs within 1 kpc 
of the black hole. For the last two, we show e.g. the gas mass 
within spheres of five different radii, equally spaced from 200 pc 
to 1000 pc. Careful comparison of the ramses and gasoline snap¬ 
shots in Figure[2and in the top panels of Figure [^indicates some 
differences in the galaxy orbits (e.g. in ramses the galaxies are 
farther apart at first apocenter), and differences in gas structures 
(e.g. in GASOLINE gas forms into smaller clouds). We will address 
these differences in §3.2| 

Overall, however, the dynamics of the black holes is similar 
until we can follow them in a consistent way. We note here that in 
our GASOLINE simulations the black holes never merge, whereas in 
RAMSES they merge when the galaxies coalesce, shortly after third 
pericenter, once they pass within each others’ accretion radii (in 
this case about 60 pc). Black hole mergers are implemented in 
GASOLINE, but we chose to turn them off to study the dynamical 
behavior ( Van Wassenhove et al.|2014t|Capelo et al.|2015| ). This 
choice affects accretion on the black holes, and therefore we re¬ 
frain from commenting on black hole growth at late times in the 
simulation. 

The patterns in star formation and black hole accretion 
driven by dynamics, i.e., how merger-driven inflows at pericen- 
ters enhance them, are also similar, once the differences in the 
orbits are taken into account. Indeed, one of the main conclu¬ 
sions of our experiment is that even when trying to make the 
set-up of our simulations as close as possible in the initial condi¬ 
tions and in the choices of implementations, the way parameters 
are set up is intrinsically different, and therefore differences arise 
even in a controlled experiment. 

Qualitatively, the behavior in gasoline and ramses simula¬ 
tions is similar. At first pericenter, the black holes pass within 
about 10 kpc of each other, but this induces little gas inflow and 
little change in the BHARs or SFRs. Between first and second 
pericenter, the galaxies act essentially as though they are iso¬ 
lated, with relatively steady SFRs and BHARs. During this time 
the SFRs and especially the BHARs fluctuate on short timescales 
- in the BHAR case, fluctuations of an order of magnitude are 
common. These fluctuations result from stochastic fueling of the 
central black hole ( [Hopkins & Hemquist|2006| ), driven both by 
structure in the interstellar medium (ISM) and by AGN feedback 
( Gabor & Boumaud|2013[ ). Following [Capelo et al.[ ( [2015] ), we 
call this the stochastic phase of the merger. 

The action begins at second pericenter. As the two black 
holes pass within a few kpc, tidal torques trigger an increase in 
gas mass within the central few hundred pc in both the primary 
and secondary galaxies. The central gas mass increase, more pro¬ 
nounced in the secondary, triggers enhanced central star forma¬ 
tion and black hole growth. At third pericenter ~ 150 Myr later. 
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Fig. 3. Face-on images of gas surface density of an isolated disk galaxy simulation, using gasoline (left) and ramses (right). The characteristic 
sizes of gas structures are larger in ramses than in gasoline. This owes partly to the Jeans polytrope pressure floor implemented in ramses, and 
partly to differences in more fundamental aspects of the codes (e.g. hydrodynamic or gravity solvers). 



Fig. 4. Power spectra of gas surface density fluctuations in isolated 
GASOLINE and ramses galaxies. The power spectra are normalized by 
a third-degree power law to highlight the differences between simu¬ 
lations. GASOLINE shows more power at small scales, and less power at 
large scales, emphasizing that gas structures in our gasoline simulations 
tend to have smaller sizes than those in ramses. 

the interaction induces another burst of activity. In both simula¬ 
tions, this picture (starting at second pericenter) roughly follows 
the well-known scenario of merger-driven starburst and AGN fu¬ 
eling ([Sanders et aL| 19881 |Barnes & Hemquist|1991[|Di Matteo 
jet al.|2005D . 

In summary, the overall behavior of the ramses and gaso¬ 
line simulations is similar, and broadly consistent with previous 
studies of galaxy mergers. There are, however, important dis¬ 
crepancies between the two simulations, which we address next. 
Understanding these differences and why they arise is important 
in order to extract and retain the results that we can consider 
robust. 
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3.2. Quantitative discrepancies 
3.2.1. Orbits 

The galaxy orbits differ slightly between the gasoline and ram¬ 
ses runs, as implied by the top panels of Figure In the gaso¬ 
line run, the galaxies move to wider separations between first 
and second pericentric passages (around t = 500 Myr), and the 
second pericenter occurs at a later time (~ 1.0 Gyr rather than 
-0.8 Gyr in the ramses case). Since the mergers are initialized 
with the same positions and velocities in both simulations, we 
attribute these differences to mass discrepancies. 

As described in ^ galaxy halos in ramses simulations must 
be truncated to fit inside a computationally reasonable box size. 
When we created the initial conditions, we chose to keep the total 
halo mass constant, so the truncation increases the dark matter 
density by a few percent inside the truncation radius. If instead 
we had kept the central density constant, then the total halo mass 
would be modified and the interaction orbit would start to dif¬ 
fer even at large separations, before the two galaxies overlap or 
develop tidal features. Instead, our choice keeps the early-stage 
configuration of the interaction relatively unchanged (see Figure 
[^. Subsequently, the higher central density of the truncated ram¬ 
ses halo increases the gravitational forces at shorter distances, 
once the two haloes start to significantly overlap. We attribute the 
faster merging in ramses to this difference, rather than to effects 
of the codes themselves (e.g., Poisson solver type or accuracy). 
To further probe this, we estimated analytically the timescale re¬ 
quired to “free-fall” from the first apocenter to the next encounter 
at < 5 kpc, assuming the radial profile of the haloes and bary- 
onic components did not evolve from the initial conditions: the 
timescales are 382 Myr for the gasoline initial mass distribution, 
and 327 Myr for the ramses one. The timescales measured in the 
simulations are about 350 and 290 Myr, respectively: they are 
both shorter than our analytic estimates (probably due to the ex¬ 
tra effects of dynamical friction, and/or increased mass concen¬ 
tration at the first pericenter), but the relative difference of ^15% 
is fully consistent with that in the analytical estimate. Hence, the 
faster merging timescale in the ramses merger is fully consistent 
with having resulted from our choice of keeping the total halo 
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GASOLINE RAMSES 



Log(Density [cm^]) Log(Densrly [cm^]) 

Fig. 5. Density-Temperature diagrams (top row) and density PDFs (bottom row) for gas in gasoline (left column) and ramses (right column). 
We show only gas within 500 pc of the primary galactic center in each simulation, and we use a snapshot near the maximum separation (first 
apocenter, around 500 Myr). In the top panels, we include the temperature distribution along the left axis (arbitrarily normalized, and shown on 
a linear scale). We also schematically show the Jeans polytrope temperature floor for the ramses simulation (straight black line), which keeps the 
average gas significantly hotter than in gasoline. 



Fig. 6. Toomre Q parameter estimated as a function of radius for the 
initial disk galaxies in a 4:1 merger in ramses. A horizontal line marks 
the marginally-stable value g = 1. At a gas temperature of 10^ K (which 
is effectively enforced by the Jeans polytrope temperature minimum), 
the secondary galaxy is substantially more stable than the primary, and 
it remains stable or near marginal stability at all radii. This leads to a 
low SFR in the secondary. In the gasoline simulations without a Jeans 
poly trope, the secondary has higher SFRs. 


mass constant when a truncation is applied to the initial condi¬ 
tions. 


3.2.2. Black hole accretion rates 

During the relatively quiescent phase between first and second 
pericenter, the BHARs differ significanly in the two simulations 
(see Figure]^: ramses BHARs fluctuate around 10“^-^ M© yr“^, 
while those in gasoline fluctuate around 10“^-^ M© yr“^. In any 
case, the black hole growth during this period is small because 
the accretion rates are low. Thus the detailed level of black hole 
fueling is relatively unimportant during the quiescent phase. 

During the merger phase, the primary BHs show marked dif¬ 
ferences in the BHAR. In gasoline, the primary BHAR peaks 
sharply (by a factor >10) just after second pericenter, remains 
elevated (with large fluctuations) more-or-less until third peri¬ 
center, and rises again just after third pericenter. In ramses, the 
primary shows a brief spike in BHAR at second pericenter, then 
returns to the quiescent level even through third pericenter, un¬ 
til experiencing another spike at coalescence. The secondary 
BHARs, in contrast, are similar in both simulations, showing a 
peak in BHAR at second pericenter and another peak at third 
pericenter that lasts through coalescence. In ramses, black hole 
coalescence occurs ~ 50 Myr after third pericenter, whereas co¬ 
alescence never occurs in gasoline. The coalescence limits dual 
black hole growth in the merger remnant. 

We attribute these quantitative discrepancies between black 
hole growth in gasoline and ramses mainly to differences in the 
black hole fueling and feedback models. In gasoline, the accre¬ 
tion rate includes a boost factor of 3 to the formal Bondi rate, 
and the feedback coupling efficiency is set to a relatively low 
value of 0.1% (calibrated from idealized mergers and zoom sim¬ 
ulations; [BeIIowyetaI]20^. Thermal energy from feedback 
is dumped into a single nearby particle at every timestep with¬ 
out any storage. In ramses, on the other hand, the accretion rate 
includes no boost factor, and the coupling efficiency is set to rel- 
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Fig. 7. Same as Figure]^ but for a 2:1 merger rather than a 4:1 merger. BH separation in kpc (top), BHAR in M© yr ^(second row), gas mass in 
M© within spheres up to 1 kpc in radius (third row), and SFR in M© yr“^ within the same spheres (bottom row) vs. time for a 2:1 merger with 
GASOLINE (left) and ramses (right). For each simulation, we show quantities for the primary galaxy/BH on the left, and the secondary galaxy/BH 
on the right. We show gas masses and SFRs within radii of 200 pc to 1 kpc, in increments of 200 pc. Vertical dotted lines mark local minima in the 
separation between the two black holes in all panels. The gasoline and ramses simulations show qualitatively similar results: enhanced accretion 
and star formation rates at second pericenter, and especially at 3rd pericenter/coalescence. 


atively high value of 0.15 (calibrated from cosmological simula¬ 
tions with much lower resolution, and slightly different feedback 
prescriptions; see [Booth & Schaye|[2009| [Dubois et al.||2012| . 
Feedback energy is dumped into the entire accretion region, but 
only if enough energy is “stored” to heat the gas in that region to 
10^ K, which increases the effective efficiency. 


Black hole growth is mostly self-regulated in these simula¬ 
tions, suggesting that the difference in coupling efficiency causes 
more of a difference than the boost factor. The efficient feed¬ 
back in RAMSES keeps the gas immediately around the black hole 
hot and diffuse, leading to lower accretion rates than in gaso¬ 
line. This effect persists throughout both the stochastic phase 
and merger phase. Owing to the efficient feedback, bursts of ac¬ 
cretion in RAMSES (see Figure are short-lived and they drive 
only minor BH growth. 

The implementation of supernova (SN) feedback in this set 
of RAMSES simulations may also contribute to suppressing black 
hole growth. [Dubois et ar] ( |2Q15[ ) find that when the SN feedback 
implementation in ramses includes delayed cooling black hole 
growth is suppressed in galaxies with bulge mass below ~ 10^ 
Mq, which is comparable to the bulge mass in our runs. 

This BH growth discrepancy raises questions about the 
choice of AGN feedback efficiency. Our best method to cali¬ 
brate the efficiency relies on comparing simulated black hole 
growth to observed BH demographics and MeH-galaxy rela¬ 
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tions, especially Mbh -^buige- This was done using large-scale 
cosmological simulations for ramses, and idealized mergers and 
zoomed simulations for gasoline ( [Bellovary et al.|2013[ ). Com¬ 
bined with the discrepancies in BH growth, this suggests that the 
optimal efficiency could depend on resolution (and stellar feed¬ 
back effects). If so, this calibration could be especially problem¬ 
atic when pushing the boundaries of simulation resolution - it is 
computationally expensive to run several high-resolution simu¬ 
lations just for calibration. 

In summary, different physical prescriptions for black hole 
fueling and feed back dominate over differences in hydrody¬ 
namic method (cf. [Hayward et al. 


2014). The efficient AGN feed¬ 


back in RAMSES leads to lower B 


i accretion rates than in gaso¬ 


line. The Mgas plots in Figure |2 - which show less gas inflow 
at late stages in ramses than in gasoline- suggest that gas dy¬ 
namics play a role as well. In the next sections, we explore how 
the models for gas thermodynamics influence the gas structure 
of the galaxies, which in turn will influence the details of gas 
dynamics. 


3.2.3. Galaxy gas structure 

Figure (especially the middle panels) suggests a key differ¬ 
ence in gas structure between our ramses and gasoline simula¬ 
tions: gas forms larger, smoother structures in ramses. In order 
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to eliminate any possible effects of the galaxy merger, we com¬ 
pare these same primary galaxies but in an isolated context. We 
show images from these in Figure 

In the case of gasoline, we show a snapshot at 200 Myr of 
the 4:1 merger. This is well before the first pericenter, and thus 
the galaxy appears as it would in isolation. For the ramses case, 
we recall that the galaxies do not undergo a separate relaxation 
simulation before beginning the merger simulation, as do the 
GASOLINE galaxies. Thus the galaxies relax during the first 200- 
300 Myr before first pericenter. In order to make a fair compari¬ 
son with the isolated gasoline galaxy, we ran a separate ramses 
simulation with an isolated version of the primary galaxy. We 
show a snapshot from this isolated ramses simulation at about 
250 Myrs, which is similar to the “age” of the gasoline galaxy. 

Figure [^reinforces the difference in gas structure. In ramses 
the gas forms thick filaments and large clumps, whereas in gaso¬ 
line the gas forms many smaller filaments and clouds, with a 
more fiocculent appearance. Moreover, dense gas structures ap¬ 
pear to extend to larger radii in gasoline, while the ramses outer 
disk is quite smooth. 

To quantify the differences in gas structure, we calculate 
the power spectrum of two-dimensional spatial variations (as in 
[Christensen et al.|2012|). We first crop each gas surface density 
map shown in Figure to a width oi L ^ 6 kpc to isolate the 
galaxy ISM, then calculate the 2D fast Fourier transform. This 
yields a 2D image in frequency space, 2(i/,v), where u and v 
are frequency coordinates. We shift the image so that zero fre¬ 
quency is at the image center, and calculate the 2D power as 
P(u,v) = Pp. From the 2D power image we calculate the ID 
power spectrum, P(k), by calculating the average power per pixel 
in a concentric set of circular annuli. Here, k refers to the radius, 
in pixels, of each annulus placed on the frequency-space image; 
k is a spatial frequency coordinate related to spatial wavelength 
A via A = LIk. Structures of size 5' correspond to a wavelength 
of A = 2s. Finally, we normalize the power spectrum by a power 
law to clarify the differences between simulations: in Figure]^ we 
show P'(s) = F(5’)/(0.0245’^). a third-degree power law roughly 
fits the shape of P(s), so this normalization flattens the trend. 

Figure]^ shows that our gasoline simulation has more power 
than RAMSES on scales 5' smaller than 100 pc, and less power 
on scales larger than ~ 200 pc. This confirms that the sizes of 
gas structures in gasoline indeed tend to be smaller than those 
in RAMSES. Both simulations have the same amount of power 
on scales ~ 100 - 200 pc, which is approximately the width of 
filaments and clumps in the ramses simulation (i.e. the Jeans 
length). The gasoline simulation includes structures of this size, 
but their internal fragmentation leads to the enhanced power on 
smaller scales. 

In the following section, we argue that the relative excess of 
gas structure on large scales in ramses is due to the implementa¬ 
tion of a Jeans polytrope, the pressure floor which stabilizes the 
gas (which is absent from the gasoline simulations). 

3.2.4. Stabilization by the Jeans polytrope 

In Figure we show density-temperature diagrams and den¬ 
sity distributions in both gasoline and ramses. We show only 
gas within 500 pc of the primary BH at a timestep during the 
stochastic phase, between first and second pericenter. This gas is 
representative of gas in the disk of the galaxy. 

We also schematically show the Jeans polytrope in the ram¬ 
ses density-temperature diagram. By construction, gas in ram¬ 
ses is not allowed to fall below this line. In contrast, a signif¬ 


icant amount of gas in gasoline exists at lower temperatures 
and higher densities than imposed by this floor in ramses. At 
the GASOLINE temperature floor of 500 K, gas at the threshold for 
star formation (100 cm“^) is resolved by 64 gas particles (with 
poorer relative resolution with increasing gas density). We note 
that GASOLINE includes an option to use a Jeans polytrope, but it 
frequently goes unused in galaxy formation simulations (as here) 
- mainly because the cold dense gas to which it applies is always 
star forming gas, which is treated with a sub-grid star formation 
model. 

The thermal pressure imposed by the effective density- 
dependent temperature floor in ramses helps stabilize the 
galaxy’s gas disk against gravitational collapse. Gas does not 
reach densities higher than about 10^ cm“^, whereas the high- 
density tail in gasoline, where the temperature floor is inde¬ 
pendent of density, reaches 10"^ cm“^ (bottom panels of Figure 
[^. The gas temperature distribution (along the left axis of the 
density-temperature diagram) in ramses has a broad peak around 
10"^ K, while in gasoline it shows a bimodal structure: a small 
peak around 10"^ K, and a larger peak around 500 K (the tem¬ 
perature floor). In both simulations a substantial portion of the 
gas reaches the temperature floor, but in ramses the temperature 
floor is set by the Jeans polytrope. 

The higher typical gas temperature in ramses implies a more 
pressurized ISM, and a larger Jeans length (and Jeans mass). The 
larger Jeans length naturally leads to larger collapsed structures, 
giving rise to the larger filaments and gas clumps seen in Fig¬ 
ures and This is reflected by the relative enhancement of 
power at scales > 200 pc in the gas power spectrum (Figure]^. 
The higher temperature in ramses also affects the overall disk 
stability, which we address in |3.2.5| 

To test whether the Jeans polytrope is the main driver of dif¬ 
ferences in gas structure, we ran ramses simulations of an iso¬ 
lated galaxy without the Jeans poly trope, and repeated the power 
spectrum analysis of |3.2.3| In this case, the power discrepancy 
on scales > 200 pc between ramses and gasoline disappears (not 
shown). This implies that the Jeans polytrope is enhancing large- 
scale power by increasing the Jeans length. The relative deficit of 
small scale power (< 200 pc) in ramses, on the other hand, per¬ 
sists even without the Jeans polytrope pressure floor (and also 
with changes to feedback models). The small-scale discrepency 
thus seems to arise from more fundamental aspects of the sim¬ 
ulations: numerical diffusion between gas cells in ramses could 
smooth out small-scale perturbations, or differences in the grav¬ 
ity solvers could lead to different levels of clumping. 

Large gaseous clumps in ramses can have masses >10 times 
the black hole mass. The largest of these clumps may also con¬ 
tribute to the black hole scattering discussed in |2.2| We have 
performed a simple analytical estimate of the scattering induced 
by the clumps on the black hole motion on one ramses output. 
We calculated the change in the black hole’s velocity both due 
to distant encounters ( jBinney & Tremain^|1987| ) and assuming 
a single impulsive scatter that conserves energy and momen¬ 
tum, thus obtaining a lower and upper limit to the dynamical 
effect. In the case we analyzed the ensemble of the clumps was 
not producing any strong effect, but a single large clump with a 
mass 200 times larger than the black hole, and very close to it 
(~ 60 pc), could, in principle, have imparted a kick sufficient to 
displace the BH from the center. 

3.2.5. SFRs 

Overall SFRs during the stochastic phase are similar in both sim¬ 
ulations. In both cases, the total SFRs are dominated by those in 
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the primary galaxies. This dominance is more pronounced in the 
RAMSES simulation, as the SFR of the secondary galaxy is quite 
low (see Figure]^. Notably, the gas mass within 1 kpc is similar 
for the two simulations, implying that ramses has a lower star 
formation efficiency (SFR/Mgas). 

We attribute the abnormally low SFR in the secondary ram¬ 
ses galaxy - but not in the primary galaxy - to the pressure sup¬ 
port provided by the Jeans polytrope. As described above, the 
Jeans polytrope acts as a minimum temperature that ensures the 
Jeans length is always resolved. 

From Figure the ramses gas density distribution peaks 
around ~ 10^ cmf^, at which density the temperature floor is 
around 10^ K. We can quantify the stability imposed using an 
analysis of the Toomre Q parameter ( |Toomre|[l964| ). In a thin 
galactic disk, the value of Q quantifies the stability of the disk 
against gravitational collapse: small values indicate gravitational 
collapse, large values indicate stability, and 2 ~ 1 indicates 
marginal instability. We estimate 2 as a function of radius us¬ 
ing: 


TiGfL 


( 1 ) 


where Vcirc/^ is the orbital frequency, Vcirc(2 is the circular ve¬ 
locity at radius r, Cs is the gas thermal sound speed, G is the grav¬ 
itational constant, and E is the gas surface density in the disk at 
rj^We calculate Cs assuming a gas temperature of 10^ K, which 
is roughly the limit imposed by the Jeans polytrope in ramses. 

Figureshows 2 as a function of radius for the initial galax¬ 
ies in the 4:1 merger. The primary galaxy is Toomre-unstable 
within a radius of ~ 2.5 kpc. The smaller, secondary galaxy is 
significantly more stable than the primary, and Q barely reaches 
below 1 (and only in the central kpc). Thus, the Jeans polytrope 
stabilizes the smaller gas disk, inhibiting gravitational collapse 
and star formation. 

This stabilization may also help explain why the SFRs near 
the secondary black hole during the merger phase (after sec¬ 
ond pericenter) are lower in ramses than those in gasoline. The 
amount of gas within 1 kpc of the secondary BH is similar in 
both simulations, but the SFR is lower in ramses (see Figure 
[^. Another contributor to the stability at smaller radii could 
be the presence of the artificially-massive black hole, which 
would boost the circular velocity and therefore Q. In the primary 
galaxy, however, we see little evidence for a suppression of the 
SF efficiency in the central regions of the ramses simulation rela¬ 
tive to the GASOLINE simulation, suggesting that the extra-massive 
black hole plays little role. Finally, efficient AGN feedback may 
help maintain the gas at higher temperatures, suppressing the SF 
efficiency. 

Near the primary BH, at and after third pericenter, the gas 
inflow to the central kpc is more significant in gasoline than in 
RAMSES. This leads to a larger gas mass and stronger burst of SF 
in GASOLINE than in ramses (see Figure [^. The origin of this dif¬ 
ference is unclear. It could result from differences in the dynam¬ 
ics around the black holes, or from the efficient AGN feedback 
in RAMSES that evacuates gas from the nuclear regions. 

Another possible driver of differences in SFR could be the 
supernova feedback recipe. In gasoline, gas heated by super¬ 
novae is not allowed to cool for a time that depends on local 
gas conditions. In ramses, the cooling delay is fixed at 20 Myr. 
Supernova feedback helps regulate star formation, so different 

^ This estimate is accurate for flat rotation curves, for which the 
epicyclic frequency is /c = ^/2Q = V2vcirc/r. It is an under-estimate of 
Q in the central kpc or so of our simulations, where the rotation is rising. 
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SFRs could arise from these different implementations. The pri¬ 
mary galaxies have similar SFRs during the stochastic phase of 
our simulations, so the different supernova feedback seems to 
have a minor effect. 

In summary, the two simulations have similar SFRs overall, 
but the the gasoline simulation shows more SF in the secondary 
galaxy and during the merger phase than does the ramses sim¬ 
ulation. We attribute these differences to stabilization due to the 
Jeans polytrope in ramses. This temperature floor increases the 
gas stability against collapse, suppressing star formation. 


3.3. A 2:1 merger 

We show the results of a 2:1 galaxy merger in Figure [ 7 ] Broadly, 
the behavior is similar to that in the 4:1 galaxy merger. The black 
holes undergo a relatively quiescent stochastic phase of growth 
lasting until second pericenter, at which point gas inflows trigger 
more activity. The star formation rates follow a similar broad pat¬ 
tern: steady star formation until second pericenter, then a boost 
triggered by gas inflows. 

In the RAMSES simulation, the BH accretion rates during the 
stochastic phase are again lower than in gasoline, owing to dif¬ 
ferences in the supernova and AGN feedback model. The ram¬ 
ses primary BH shows little enhancement of accretion at second 
pericenter, while the gasoline primary shows a burst just after 
second pericenter. At third pericenter, both simulations show a 
burst of accretion. The secondary BHs show the opposite trend: 
the GASOLINE secondary shows little enhancement at second peri¬ 
center, and the ramses secondary shows a brief burst at second 
pericenter. 

The star formation rates are broadly consistent between the 
simulations. The secondary galaxy in this 2:1 merger is suffi¬ 
ciently unstable (owing to higher gas surface densities) in ram¬ 
ses to have comparable SFR during the stochastic phase to that 
in GASOLINE. Both simulations show SFR enhancements between 
first and second pericenter, but these occur at different times and 
locations in the two simulations. The gasoline primary shows 
an SFR enhancement about 200 Myr before second pericenter; 
the GASOLINE secondary shows an unusual SFR peak at around 
1 kpc within 100 Myrs after first pericenter, then a weak cen¬ 
tral burst about 200 Myrs before second pericenter. In ramses, 
both the primary and secondary galaxies show SFR bursts about 
100 Myr after first pericenter, but the secondary one is much 
stronger. SFR enhancements during the merger phase again ap¬ 
pear stronger in the gasoline simulation, but the differences are 
less pronounced than in the 4:1 simulation. 


4. Conclusion 


We have compared galaxy merger simulations including black 
hole growth with ramses, an AMR code, and gasoline, an SPH 
code, to find what is robust against differences in the codes 
and their physics implementations. We ran a small suite of new 
RAMSES simulations and used gasoline simulations from a suite 
described in [Capelo et al.| ( [2015| ). In both cases, the simula¬ 
tions had a spatial resolution of ~ 20 pc and mass resolution of 
~ 5 X 10^ Mq, and included standard sub-grid recipes for gas 
cooling, star formation, supernova feedback, and black hole fu¬ 
eling and feedback. 

The simulations show general agreement in their dynami¬ 
cal behavior, and the results are broadly in line with the long¬ 
standing picture of galaxy merger simulation evolution (cf. [PI 
[Matteo et al.||2005[ [Hopldns et al.||2006| ). Black hole growth is 
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relatively slow while the galaxies are well-separated, including 
the time between first and second pericenter. During this phase, 
black hole growth is driven by stochastic a ccretion ([Hopkins &| 
|Hernquist||2006[ [Gabor & Bournaud|[2013[|Capelo et al.||2015| ), 

just as for isolated galaxies. Once the galaxies’ effective radii 
overlap at second pericenter, tidal torques trigger gas inflows 
toward the nucleii. The activity is strongest when the galaxies 
finally coalesce (beginning at third pericenter). The inflowing 
gas, driven to higher densities, initiates peaks in the SFR and the 
black hole accretion rate. 

While the general behavior of the galaxy mergers is similar, 
there are quantitative differences between simulations run with 
different codes. We attribute most of these discrepancies to dif¬ 
ferences in sub-grid cooling and feedback models, but effects 
of the hydrodynamic methods (SPH and AMR) or even grav¬ 
ity solvers (gravity tree and Particle Mesh) may play a role. In 
RAMSES simulations, nuclear black holes accrete significantly less 
than in gasoline, both while the galaxies are well-separated and 
during coalescence. This arises due to a higher feedback effi¬ 
ciency in RAMSES, which regulates the black hole growth at a 
lower level. Furthermore, the more efficient black hole feedback 
ensures that central star formation is lower in ramses during the 
merger and coalescence. Another difference between codes is the 
ISM gas structure: in gasoline gas forms small filaments and col¬ 
lapses into tiny dense knots, whereas in ramses gas forms thick 
filaments and large clumps. Gas structure differences owe partly 
to the inclusion of a Jeans poly trope pressure floor in ramses, 
though at the smallest scales they may be related to numerical 
methods. 

While differences in hydrodynamic method are important in 
some regimes of galaxy evolution models (e.g. [Hayward et al. 


the treatment of baryonic physics (cf. [Scannapieco et al.|2012| ). 
Future work will better address the effects of hydrodynamic 
method by making the baryonic physics treatments as similar 
as possible, but due to differences in code structure, this is not 
always possible. 

Our study underscores that simulation codes that success¬ 
fully reproduce a range of observables are robust at predicting 
the same general physical behaviour in galaxy mergers, how¬ 
ever, one should be careful when comparing quantitatively one 
particular simulation to one given observable, as these can be 
model-dependent. The trends, on the other hand, are robust, and 
can be used to obtain physical insight. 
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